Self-similar solutions of the G-equation - analytic description of the flame surface 
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The main feature of the flame kinematics can be desribed with the G-equation. We investigate 
the solutions of the G-equation with the well-known self-similar Ansatz. The results are discussed 
and the method how to get self-similar solutions is briefly mentioned. 

To understand the dynamics of flame fronts is a general interest for combusion engineering. One way to gain such 
knowledge is the investigation of the G-equation model. Since the G-equation was introduced by Markstein in 1964 
■ [H numerous analytical and numerical approaches 3-5] have been used to find the solutions. The study of @ uses 
the method of characteristics for nonlinear first order partial-differential equations(PDA). 

Dekena Q uses a computaional fluid dynamics(CFD) code Fire to investigate turbulent flame propagation in 
U ■ Gasoline Direct Injection engines. Further references for application of the G-equation is given therein. 

| In the G-equation model the flame itself is treated as a surface (flame front) that separates the burnt from the 
{SJ unburnt gas. A detailed derivation of the model can be found in |2j. The flame front is described here by the relation 
1 of 



^ 9C(M) + Sjl + ( d -^) =0. (2) 



in 
d 



G(W) = *-CM=0. (1) 
The most general form of the G-equation is the following: 




c/3 ' dt dr V \ dr 
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*55 ■ Where u(r, z,t) and v(r,z,t) are the radial and the axial components of the gas velocity, and Sl is the constant 
^-y modulus of the laminar burning velocity, respectively. 

[ In the following we use a well-known method, which is very popular to investigate nonlinear PDAs, namely searching 
for self-similar solutions. We are looking for solution of ([2]) in the form of 

CM) :=t- a f(v). (3) 



The similarity exponents a and ft are of primary physical importance since a represents the rate of decay of the 
magnitude ((x, t), while ft is the rate of spread (or contraction if ft < ) of the space distribution as time goes on. 
This kind of solutions have a broad range of applicability, e.g. a large number of diffusion- reaction- type equations @ 
or time-dependent telegraph-type heat-conduction equation [7j can be examined as well. 
Substituting this into ([2j we have 

- at-*- 1 /^) - ftt- a -P- x f{ri)r - u(r, z, t)^-^'^) - v(r, z, t) + (4) 



where prime denotes differentiation with respect to r\. 

After setting set u = Si = 1 and v = 0, one can see that this is a non-linear ordinary differential equation(ODE) if 
and only if a = — 1 and ft = 1 (the universality relation). The corresponding ODE we shall deal with is 

m - Vf'iv) - f'(v) + V~l + WW\ = 0- (5) 
The solutions can be easily obtained and read: 



=cir)-ci - yi + c?; f(rj) 2 = ciVw? - 2) (6) 

where c\ is the real integration constant. The solutions can be seen in Fig 1. The corresponding self-similar solutions 
are 



'r(r - 2b) 



C(r,t) 1 = cir-ci- A /l + <5i C{r,t)2 = c 1 tJ( v ^ 1 ). (7) 
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Both solutions can be prooven with direct derivation. The first solution is trivial however the non-trivial second one 
is presented on Fig. 2. 

We may consider more sophysticated flow systems where thew radaial and axial componenst of the gas velocity and 
the Sl are functions of time and radial position. First let's consider that the velocities and the constant modulus of 
the laminar burning velocity are real numbers (u, v, Si e K.) Now the ODE has the following form 



f(v) - vf'iv) - "/'fa) - « + Sty/1 + [fin) 2 ] = 0. (8) 

The solutions became a bit more complicated 

f(r))i =cir) + uc 1 +v-Si\Jl+~4; /(r?) 2 = a\J ' ~v? - rj 2 + Sl - 2r,u + v. (9) 

The function f(rf)2 is presented in Fig. 3. with the following fixed set of parameters (u = 2, v — — 1, c\ = 1, Sl = 3). 
Note, that the solution has now a compact support with a non-vanishing first derivatives at the border. This means 
that there is no flux-conservation at the boundaries. There are more analytical solutions for special u(x, t), v(x, t) and 
Si(x, t) functions. The only remaining task is to give reasonable physical interpretation for u(x, t), v(x, t) and Si(x, t). 

Another generaly interesting question is the dispersion relation and the attenuation distances. It can be examined 
how wave equations or other nonlinear evolutionary equations propagate plain waves in time and in space. Inserting 
the usual plain wave approximation C( r )*) = e l ( kr+ult ' into © the dispersion relation and the attenuation distance 
can be obtained. These are the followings 

(10) 
(11) 



Re(k) - l + k 2 \(l - signum[-l + fc 2 ] 
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Im(k) bJ + k + \y/\-l + k 2 \(l- signum[-l + fc 2 ]) 

inserting the relation of u) = kc where c is the propagation velocity of the signal the formulas only depend on the 
wavenumber vector k. Considering c = 1 propagation speed Fig. 4. shows the phase velocity as the function of the 
wavenumber. The (1 — signum[— 1 -I- k 2 ]) in the formula is responsible for the compact support of the function. Fig. 
5. shows the attenuation distance of the various waves. At fc = 1 the a(k) function is non-analytic. 
In this short study we just wanted to present that self-similar solution can be easily used to generate analytic 
solutions for the G-equation. With more general, and more physical relations for the radial and axial gas flow velocity 
hopefully more physical solutions can be obtained. Anyhow, any kind of analytical solution of a nonlinear PDA can 
be usefull giving a solid basis for testing complex and sophisticated two or three dimensional numerical finite-element 
computational fluid dynamics codes. 
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FIG. 1: Solutions of Eq. (5), red line presents f{rj)2 and the green one shows f(rj)i . 




FIG. 2: The self-similar solution C(r, t)2, from Eq. (7) in the range r = .. 15, t = -2 .. 7 for ci = 1. 
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FIG. 3: The solution 7(77)2, U = 2,v = -l,ci = 1, S L = 3. 




FIG. 4: The dispersion relation v p (k) for Eq. (2). 
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FIG. 5: The attenuation function a(k) for Eq. (2). 



